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INTRODUCTION 


The Grant NGR-06-002-1A7 , "Theoretical Investigations on Plasma 
Processes," is concerned with i) the system analysis of plasma centrifuges 
for isotope separation, and ii) the deposition of the sputtered atoms on 
system components such as the solar energy collectors of ion propulsion 
systems. The progress made on these subjects in the period from 9.1.75 to 
11.1.76 is communicated herein. 

The present report consists of three papers, two of which have been 
accepted for publication. By means of system analyses, it is shown that 
plasma centrifuges at high plasma and power densities are technically 
promising for isotope separation. The largest spatial isotope separation 
is obtained by maximizing the Hartmann number and the Hall coefficient of 
the isotope mixture plasma of the centrifuge discharge. Furthermore, it 
is shown that induced magnetic fields have no effect for small Hall- 
coefficients and little effect for large Hall-coefficients on the plasma 
rotation and the isotope separation efficiency. A new approach to the 
deposition problem of sputtered atoms is given. Its solution leads to a 
simpler, but still mathematically difficult integral equation. 
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BOUNDARY-VALUE PROBLEM FOR PLASMA CENTRIFUGE 
AT ARBITRARY MAGNETIC REYNOLDS NUMBERS*) 

H. E. Wilhelm and S. H. Hong 

Department of Electrical Engineering 
Colorado State University 
Fort Collins, CO 80523 

ABSTRACT 

The boundary-value problem for the partial differential equations, 
which describe the (azimuthal) rotation velocity and induced magnetic 
fields in a cylindrical plasma centrifuge with ring electrodes of dif- 
ferent radii and an external, axial magnetic field, is solved in closed 
form. The electric field, current density, and velocity distributions 
are discussed in terms of the Hartmann number H and the magnetic Reynolds 
number R. For small Hall-coefficients, (jjt « 1, the induced magnetic 
field does not affect the plasma rotation. As a result of the Lorentz 
forces, the plasma rotates with speeds as high as 10 cm/sec around its 
axis of symmetry at typical conditions, so that the lighter (heavier) 
ion and atom components are enriched at (off) the center of the discharge 
cylinder. 

*) Supported in part by NASA. 
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On principle, electrofliagnstlc forces allow to rotate plasmas up 
to relativistic speeds. Theta pinch experiments show that the plasma 
rotates during the discharge pulse at such high speeds that the energy 
distribution of the emitted neutrons Is shifted.— From the theoreti- 
cal point of view, the basic mechanism for plasma rotation by means 
of crossed electric and magnetic fields and Lorentz forces In rare- 
fied-—— and dens e— • plasmas is understood qualitatively. Experi- 

mental evidence on Isotope separation In rotating plasmas has been 

10-11 20 21 
reported among others by Bonnevler, Gullloud, — Heller and Simon, — 

22 23-24 

James and Simpson, — and Ban and Sekiguchi Exact solutions for 

the boundary-value problems describing plasma centrifuge systems are 

not known, neither for collision-dominated nor for collisionless — — 

plasmas. 

A simple model for an electrical discharge centrifuge with an 

axial, external magnetic field is shown in Fig. 1. This plasma 

centrifuge employs electrodes of different radii R, and R (R >>R ) in 

+ — + — 

the end plates z = ±c of an electrically isolating discharge chamber of 

radius R^ so that the field lines of the current density ^ and of the 

external axial magnetic field B^ cross under a nonvanishing angle (except 

at the chamber axis). The resultant Lorentz force i x B rotates the 

0 

discharge around its axis of symmetry. In steady state, the magnetic 
body forces in the azimuthal direction are balanced by the viscous forces 
(boundary layers at the chamber walls) . As opposed to the centrifuge 
with radial electric current flow between an inner and outer cylinder 
electrode, the centrifuge scheme in Fig. 1 avoids the boundary layer 
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and losses at the inner cylinder surface. In the following, the boundary- 
value problem for this centrifuge is solved in closed form based on the 
magnetogasdynamic equations for dense isotope plasmas with negligible 
Hall effect. 
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THEORETICAL FORMULATION 

For a purely azimuthal flow field, v = {0, v(r,z),0}, the plasma 
behaves incompressible, V • v != 0. From the continuity equation, 

V • (pv) = V • Vp = 0, it follows then that the density gradient Vp is 
everywhere perpendicular to the flow field v. In accordance with the 
magnetogasdynamic equations,— Maxwells equations and Ohms law^ [which 
imply * -p^aVx(vxB) ] , the plasma in the centrifuge with homogeneous 
magnetic field (Fig. 1) is described by tbe l)Oundary-vaIue problem 
for the azimuthal velocity V„(r,z) and azimuthal induction B (r,z) fields 

U u 

(r,0,z = cylindrical coordinates) : 


r 3r 3r^ 2 

r 


2 

3 V, 


3z 


B 3B- 
o 6 


yp 3z 
o 


1 

r 3r 


3r ^ 2 , 2 

r oz 


3V, 


= - y oB 


O O oz 


where 


Vg(r,z)^^j^ = 0, - c ^ z £ + c 

o 


V.(r,z) , = 0, 

0 ■ Z=» ±C 


0 < r < R 


and 




y I 

o 

2irR 

o 




-C < z < + c 


y I 6 (r-R^) 
o X 

2tt r 


0 < r < R 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


Eqs. (1) and (2) are the azimuthal components of the equations of 
plasma motion and magnetic induction, respectively (y - viscosity, 
y^ = Tnagnetlc permeability, cJ = electrical conductivity)* The 
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0 


t 


boundary conditions (3) and (4) consider that the plasma does not slip 

at the walls r = R and z = ±c. The boundary conditions (5) and (6) 
o 

follow from Maxwell's equations for a total discharge current of jl| 
amps flowing from the ring "anode" (r=»R_j^) to the ring "cathode" (r=R_) 
of vanishing radial width [6(r-R^) = Dirac function] if I < 0 (Fig. 1). 

According to Eq. (6), the net current sustaining the discharge is 

R R 

o o 

2'Jt/ j (r,z=±c) rdr * l/ 6(r-R )dr ■= I. The pressure distribution 
0 ^ 0 " 

P“p(r,z) in the rotating plasma is determined by the r- and z-components 



(7-8) 


(9-10) 


in accordance with Maxwell's equations. In absence of the Hall effect, 

WT<<1, it is Vy^ = y (i , 0, i ). Hence, B = 0 and B = 0 because of 

0 r z r z 

the homogeneity of the boundary conditions for B^ and B^, whereas 
Bg 0 since ^ ^ 0[Bg(r,z = ±c) = (y^l/2iTr)H(r-R^) ] . Since the 
Induced magnetic field is azimuthal, B - {0, Bg, B^}, the induced 
electric field is independent of Bg, i^e* V^B - VgB^e^. 




NONDIMENSIONAL BOUNDARY-VALUE PROBLEM 


The characteriatic nondlmenslonal parameters of the raagnetogas- 
dynamic centrifuge under consideration are obtained by introducing the 
nondimen sional independent and dependent variables, 

p = r/R^ » 0 1 P 1 1 , (11) 

r, = z/c , -1 1 C 1 + 1 , (12) 

and 

V(p,C) = V^(r,z)/V^, B(p,5) = B^r,z)/B^ , (13) 

where 


V = I/2ttR B ac, B = B . (14) 

o o o o o ^ ^ 

In terms of the nondlmenslonal space variables and fields, the 

boundary-value problem defined in Eqs. (1) - (6) assumes for V(p,C) 

and B(p,c) the form: 


1 9 . 3V. V . „ 

p 9p ^P9p^ " 2 ^ 

p 

2 2 
-2 9 V H 9B 

3^2 = - R 95 

> 

(15) 

1 9 . 9B. B . - 

p 9p ^^9p^" 2 ^ 

p p 

2 9^B _ R 9V 

3^2 ^2 9; 

> 

. (16) 

where 

V(p,?)p^^ =0, 

-1 ^ + 1 

> 

(17) 

V(p,e)^^^l = 0, 

0 < p < 1 

> 

(18) 

and 


-1 ? <_ + 1 

5 

(19) 


i ^ [pB(p,5)]^^^^ = R 6(p-p^)/p, (j < p ^ 1. (20) 

with ^ 


- 7 - 


7 


H 5 (a/vi)^B^R^, N = c/R^, R h u^l/2TrR^B^ = y^oV^c > o. (21) 

The llurtinunn number H, N, nnd the magnetic Reynolds number R characterize 
the ratio o£ Lorentz to viscous forces, the geometry of the centrifuge, 
and the intensity ratio of the induced and external magnetic fields, respectively 
The linear statement, 

B(p,C) = Rp + T(p,C) , (22) 

reduces the Eqs. (16), (19), and (20) for B(p,?) to equations with 
a homogeneous boundary condition (2^) for y^(p,5); 


1 3 , 3T. T . R 


(23) 


where 

'*'(P»^)p=l =0, -1 < Cl + 1 

1 3 5(P-P^) 

^ 3^ [p'i'(p.c)l^==ii " Rl— y 2], 0 < p < 1 


(24) 

(25) 


In view of Bessel's differential equation, Z" + p Z ’ + 

m m 

2 —2 2 

(k - p m )Z =0, for cylinder functions Z (k p), partial solutions 
n m mn 

of the coupled inhomogeneous Eqs. (15) and (23) are sought in the form, 


V^(p,c) = J.(k P) £JK) 
n ± n n 


(26) 

T (p,C) = J.(k P) g„(C) 
n inn 


(27) 


where the eigen-values k > 0 are determined by the homogeneous 

n 

boundary conditions (17) and (24) as the real roots of the transcen- 
dental equation. 


J, (k ) = 0 , n = 1,2,3,.... 
I n 


(28) 


Thus, the general solution of the coupled Eqs. (15) and (23) obtains 

OQ 

by linear superposition as the Fourier-Bessel series:— 
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V(p,C) = Z J,(k P) 
n-1 ^ 

00 

f(p,c) ” 2 

n=l 


(29) 


(30) 


Substitution of Eqs. (26) - (27) Into Eqs. (15) and (23) yields 
ordinary coupled differential equations of second order for f^(?) and 


f " 

n 

n 


f = 
n 

-hVr"^ g' 

n 

f 

(31) 


n 


Sn = 

-Rf ' 
n 

» 

(32) 


By elimination, Eqs. (31) - (32) are reduced to decoupled differentail 
equations of fourth order, 


f " (2k^ + H^)N"^f ‘ f = 0 

n n n n n 

g''”- (2k^ + H^)N^g;' + k^* N‘ g^ = 0 
n n n n n 

with 


(33) 

(34) 


^n<'>C=±l = ° - <35) 

8 „< 5>?.±1 = • 06 ) 

as boundary conditions, by Eqs. (18) and (25), respectively. In 

deriving Eq. (36), the Dirac function in Eq. (25) has been expanded 

28 

in the Fourier-Dini series, 


6(P - Pj/P=2 + 2 E [J (k p,)/jJ(k^)]J^(k p) . (37) 

± ,on±onon 

n=l 

In addition to Eqs,. (35) - (36), f (c) and g (<;) have to satisfy also 

n n 

the coupled Eqs. (31) - (32). With the four real roots of Eqs. (33) - 
(34) [f^, g^ exp(MC)], 
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w - U ) , CO 

In n 2n 


= CO , CO„ = -CO , CO, = -CO 

n ’ 3n n 4n n 


where 


(of 2 2"'^ N{(2k^ + H^)±[(2k^ + - 4k;;r} 


2 v 2 ,,A.k^h 

h 


( 38 ) 


(39) 


the general solutions for f^(c) and g^(C) of Eqs. (33) - (34) can be 
written as; 


n n 


+ A 
n 


n n 


+ G 

n 


sinh 


slnh 

+ 

CO 

n 

sinh 

CO K 

n ' 

sinh 

CO 

n 

sinh 

+ 

sinh 

+ 

CO 

n 

sinh 

CO C 

sinh 

‘^n 


+ B 


+ B 


n 


+ D 


+ D 


n 


cosh 


cosh 

+ 

CO 

n 

cosh 

V 

cosh 

CO 

n 

cosh 

+ 

cosh 

+ 

CO 

n 

cosh 

n 

cosh 

CO 

n 


(40) 


(41) 


Only four of the eight integration constants A^, 


are independent. 


Substitution of Eqs. (40) - (41) into Eq. (31) and Eq. (32) yields 

+ +. 2 2 2 ,+ 2 2 -1 + + 

A" [(or) - k N 1/co" = -H N R tghco" ’ D“ 

n*- n n^n nn 

B" [(to-)^ - kV]/co" = cth CO" • C" 

n *• n n n n n 


and 


+ + 2 2 2 ,+ ± ± 

C" [(co") - kV]/co“ = - R tgh co^ • B^ 

n n n n ti n 


+ +2 22 ,+ i i 

D" [(co") - k N ]/co" = - R cth CO - A 

n n n n n 


(42) 

(43) 

(44) 

(45) 


respectively. The coefficient determinant of Eqs. (42) and (45) or 

Eqs. (43) and (44) vanishes (condition for existence of nontrivial solution), 

(46) 


,+ r / v 2 „ 2 t 2 „ 2.,2 , ±.2 ^ 

A = [(co) -kNl ~HN(co) -0 

' n' n n' 


EEPKODUCIBILITY OF THE 
0EK3CNAL PAGE IS 


HO- 



10 


in agreement with Eq. (39). From the latter or Eq. (46) one deduces 


the relations. 


[(w")^ - kV]/w" «= ±NH 
n n n 


( 47 ) 


which simplify the left sides of Eqs. (42) - (45). 

Application of the boundary conditions (35) to Eq. (40) shows that 


-A = +A = A 


n 


n n 


- B” - S 
n n o 


(48) 


Substitution of Eq. (48) into Eqs. (40) and (41) gives 

4- 

sinh 0 ) C sinh w C 

f (r) = A [ 5 —] 

n'’ n , , + - 


sinh 0 ) 


n 


sinh ti) 


n 


(49) 


cosh 0 ) C cosh 0 ) C 

+ B_,[— 


cosh 0 ) 


n 


cosh u 


n 


and 


cosh, 0 ) c cosh oj C 

= -*n H t— ^ 


sinh w 


n 


sinh (D 


n 


(50) 


„ sinh w 5 sinh w 4 

s i— f - ^ 


cosh w 


cosh 03 


n n 

the latter under consideration of Eqs, (44) - (45) and Eq. (47). 
Application of the boundary conditions (36) to Eq. (50) yields, upon 
elimination, 

J^(k p ) + X (Ic p ) 

o n— on + 


A = 
n 


ra 

^n (cth m"*" + cth w ) J^(k ) 
n non 


B 

n k 


J (k p ) - J (k p , ) 
o n'^- o n + 


n (tghtj + tgh 0 )^) Jq(\) 

By combining Eqs. (49) - (52), we obtain the solutions for 

f (C) and g (C) in final form: 
n n 


(51) 


(52) 
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J (k p ) + J^(k^pJ 
r /w \ /%TTT o n o n I 

f CC)/NH T ^ . . 

(cth 03 + cth 03 )k J (k ) slnh 03 


sinh 03 C sirih 03 C 

i ~ — i-i 


n 


n n o n 


n 


slnh 03 


on- o n T 


+ - 2 + 
(tgh 03 + tgh 03 )k J (k ) cosh 03 
' “ n ° n' n o n n 


cosh 03 i; cosh 03 ^ 

[ ^ , (53) 


cosh 03 


n 


and 


J (k^p ) + J (k P.) 

/ /n I on*" on* 

g (?)/R = + V T ^ + 

(cth 03 + cth 03 )k J (k ) slnh o3 


cosh 03 ? cosh 03 Z 

J — -4- + - 


n 


n n o ' n 


n 


slnh 03 


n 


J (k p ) - J (k p , ) 
o n^- o n + 


+ - 2 + 

(tgh 03 + tgh 05 )k J (k ) cosh 03 

n n' n o n n 


slnh 03 C slnh 03 C 
[ + 


cosh 03 


n 


(54) 


Below, also the i; -derivative of is required, which Is given by 


J (k^p_) + 

g'(C)/R = + 

(cth 03 + cth u3 )k J (k ) 


"i" — 

, slnh 03 C slnh o) ^ 

r + n . - n’, 

[jjj 4- m j 

n . , + n - 


n 


n n o n 


slnh 03 


n 


slnh 05 


n 


(55) 


(tgh + tgh u‘)k^ J^(k^) 


, cosh 03 C cosh 03 C 

[y 4- (jj — J 

n , + n , - 

cosh 0 ) cosh 03 

n n 


In terms of f (?), g (O, and g’(?), the solutions for the non- 


n 


n 


n 


dimensional fields ^ = {0,V,0}, B = t0,B,l}, J - Cf^,0,J^} and E ~ 
{Ep,0,E^} of the plasma centrifuge are by Eqs. (22), and (29) - (30): 


V(p,0 = 2 J^(k^p)f^(0 

n=l 


(56) 


B(p,C) = Rp + S 
n=l 


(57) 


J (p,0 = - R"V^ S 

P n=l 


(58) 


J (P,0 = 2 + r"^ e ^ 

^ n=l 


(59) 
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The reference values V and B for V(p,C) and B(p,C) are defined in 

o o 

Eq. (14). The non-dimensional fields J (p,?) and E (p,Z) are 

P j <a P » 

normalized with respect to 


j = I/2 itR^, E 5 V B = l/2-nR oc . (61) 

o o o o o o 

If the cathode is in the plane z = -c(C = ~1) and the anode is 
in the plane z = +c(C = +1), then the reference fields and 

[Eqs. (14), (61)] are negative, since I < 0. The results are also 
applicable to the case where the anode is in the plane z = -c(C = “1) 
and the cathode is in the plane z = +c(5 = +1) ■ In the latter situ- 
ation, the reference fields V ,j »and E [Eqs. (14), (61)] are positive, 

' o o o 

since I > 0. These explanations hold for magnetic fields pointing in 
the positive z- direction, > 0; changes its sign with the sign 
of B [Eq. (14)]. Note that the magnetic Reynolds number R in Eq. (21) 
is defined to change its sign with the sign of V^. 
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AVPLTCATIONS 

As an illustration, the radial (p) dependence of the nonditnensional 
discharge fields V(p,c), B(p,^), E (p,C), J (p,?), and J (p,C) has 

M P '9 

been computed for I < 0 in the cross-sectional planes ^ = -0,99 

(cathode region) , C = 0 (central region), and C = +0.99 (anode region) 

based on Eqs. (56)- (60). The remaining field E^(p,?) is proportional 

to J (p,?) [Eq. (60)]. The characteristic (nondimensional) magnetic 

interaction number H is treated as a parameter: H = 1, 10, 100. The 

geometry parameter N = c/R is taken to be N = 1 corresponding to 

0 

= c [Eq. (21)]. The radial positions of the cathode and anode are 
assumed to be: 

p =0.01 (R = 0.01 R ) ; p^ = 0.9 (R, = 0.9 R ). 

— — o + + o 

With the exception of = B B, the dimensional fields are negative 

0 o 

everywhere where the nondimensional fields are positive, and vice-versa 

since V < 0, j <0 and E <0 for I < 0 [Eqs. (14), (61)]. 
000 

The Eqs. (56)- (60) indicate that the velocity field V(p,?), the 
current density field ^(p,C), and the electric field E^ ^(p,C) ate 

independent of the magnetic Reynolds number R, whereas the induced 

magnetic field B(p,C) is proportional to R. This is due to the azimuthal 

direction of the induced magnetic field B(p,0, which is parallel to the 

velocity field V(p,C) of rotation [Eqs . (9) - (10) ] . Accordingly, the plasma 

fields V(p,c), B(p,c)/R, J (p,?), and E_(p,?) depend only on the 

P ) ^ P 5 b 

Hartmann number H, presumed that the Hall effect is negligible (wsje B[/m 
and T are the gyration frequency and collision frequency of the electrons), 
0)T « 1 • 

- 14 - 
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Cc nitra] Region, C =* 0. In Figs. 2-6, V(p,0), [li(p,0) - Rpl/R, 

E (p,0), J (p,0), and J (p,0) « E (p,0) are shown versus 0 £ p £ 1 
P P b S 

with H = 1, 10, 100 as a parameter. It Is seen that jvj Increases 

considerably at any point 0 < p < 1 as H is increased. Similarly, 

(B - Rp)/R and the sources J ^ of the magnetic induction increase 

in intensity within the main central region 0 < p < 1 - Ap as H 

i.s Lncreased. For large values H > 10, B and J decrease in the 

P , ^ 

wall region Ap = Ap(H), so that the electrical discharge becomes more 
concentrated in the center 0 < p < 1 -Ap of the centrifuge. The 
intensity of E^ Increases uniformly in the region 0 < p 1 as H is 
increased, while 

Cathode Region, g = -0.99; The Figs. 7-11 show V(p,-0.99), 
[B(p,-0.99) - Rp]/R, Ep(p,-0.99), J^Cp, -0.99), and J^(p,-0.99) « 
Ej^(p,-0.99) versus 0 £ p £ 1 for H = 1, 10, 100. The fields V, 
and Jp ^ increase in Intensity at any point 0 < p < 1 with increasing 
H, whereas B/R decreases in 0 < p < 1 with increasing H. Since the 
ring cathode is at p = 0.01 (? = -1), the field distributions are 
closer concentrated at the axis p = 0 than those in the plane ? = 0 
(Figs. 2-6). Note that the plasma rotates only in the region p = 0.1 
with a significant velocity, since the Lorentz force decreases 

rapidly with increasing p ■+ 1. 

Anode Region, g = +0.99: The Figs. 12-16 present V(p,+0.99), 

[B(p,+0.99) - Rpl/R, Ep(p,-f0.99), Jp(p,+0.99), J^(p,+0.99) E^(p,+o.99) 

versus 0 £ p £ 1 for H = 1, 10, 100. The velocity field is fully 
developed nearly through the entire centrifuge across section 0 < p < 0.9, 
since the Lorentz force is strongest in the vicinity p = 0.9 of 

the ring anode p = 0.9(? = +1) . As a result , a thin and steep boundary 
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layer exists close to the cylinder wall (p = 1) with pJastiia counter- 
rotation at sufficiently small H-values, The radial distributions of 

B, E , J clearly indicate that, in the plane r, - +0.99, the electrical 
P ♦ C P > ? 

discharge has shifted to the region p = 0.9 due to the influence of the 

(nearby) ring anode at p = 0.9(? = +1). 

In the graphical illustrations, the cathode radius R was chosen 

to be small compared to the anode radius R^ to ensure a large angle 

between the current field lines j (r) and the external magnetic field B , 

o 

i.e. a significant Lorentz force. A comparison of the Figs. 2 and 7 

with Fig. 12 indicates that this choice of electrode radii results in 

n radial boundary layer of large width and low velocity in the lower 

half - c < z < 0 of the centrifuge. Hence, R << R (or R >> R ) is 

not the best choice for a centrifuge of maximum efficiency. Fig. 12 

demonstrates that a velocity profile raising uniformly with radius r 

and decreasing rapidly in a steep boundary layer of narrow width Ar, 

is obtained by using a cathode and an anode of the same radius R = R. < R 

- + ~ o 

which is nearly as large as the centrifuge radius R^. Although it is 
R - R , in this case, the current field lines T(r) intersect with B 
under a sufficiently large angle 0 due to the repulsion of 

the current filaments. As a result, a net Lorentz torque results for 
a centrifuge with R_ = R^ which is still of the same order of magnitude 
as for a centrifuge with R_ « R_^ (presumed that I, and B^, c, and R^ 
are the same) . 

The accuracy of the Figs. 2-16 is determined by the number of 
terms considered in the Fourier series on the computer and the accuracy 
of the eigenvalues k^. The Fourier series solutions were summed 
numerically up to n =100, and the eigenvalues k^, n = 1, 2,3 , . . . 100 , 
were computed up to the 10th decimal point. 
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The centrifuge analysis presented Indicates that extremely high 
speeds of plasma rotatiOTi are obtainable already at moderate discharge 
currents I and magnetic inductions presumed the Hartmann number H 
is not small, H > 1. As an example, consider a centrifuge discharge 
with: 

jl| = 10^ amp , |b^| = 10^ Tesla, 

a » 10^ mho m R = c = 10 ^m. 
o 

Hence, by Eq. (14) 

V = I/2itR B ac = (5 /tt) x 10^ m sec 
o o o 

and, by Fig. 2, 

0[V„] = 0[V V] = 10^ m sec”^, for H = 100. 

0 o 

If the working gas of the centrifuge discharge consists of two isotope 
gases, then the centrifugal forces would concentrate the lighter isotope 
ions and atoms in the central region and enrich the heavier isotope atoms 
and ions in the peripheral region of the discharge. According to the 
equations of motion for two isotopes of masses m, and m,, the isotope 
density ratio at distances 0 < r < -4r, where Ar is the viscous 

boundary layer thickness, is approximately (T^ = temperature of the 
isotope ions) 


n^(r) ^ n^(0) +igAm^^v(r)^/kT^ 
n^(r)'n.(0)' 




where the bar designates a spatial average over the region |z| < c. As 

a specific example, consider an uranium plasma centrifuge containing the 

O'M 1 

isotope ions (i) U and (j ) U at a temperature T^ = 10 °K (and 

electrons at a temperature T > T ) . In this case, one has Am.. = 

e o xj 

-27 -20 

m(237) - m(235) = 3.320 x 10- 'kg, kT = 1.381 x 10 ^ Joule. Hence, 
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the isotope separation ratio is: 


^ n — —1 

_£££ / _£i£ . 1.128 x 10^ for v(r) = 1 x lO^ m sec 


* 1.617 X 10^ for v(r) = 2 x lo^ in sec ^ , 

“ 2.950 X 10^ for v(r) = 3 x 10^ m sec ^ 

Based on these examples, one can assume with some confidence that high- 

power plasma centrifuges are technically realizable employing dense, 

collision-dominated isotope plasmas. The separation of isotopes by 

centrifugal forces in low density plasmas has been established experi- 
lQ-11, 20-24. 


mentally. 
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PLASM ROTATION BY LORENTZ FORGES 
IN AN ELECTRICAL DISCHARGE CENTRIFUGE WITH HALL EFFECT* 

H. E. Wilhelm and S. H. Hong 

Department of Electrical Engineering 
Colorado State University 
Fort Collins, Colorado 80523 

A system analysis for an electrical discharge centrifuge is presented 
in which the plasma rotates under the influence of the Lorentz forces 
due to the interaction of the nonuniform current density fields with an 
axial external magnetic field. The associated boundary-value problem 
for the coupled partial differential equations, which describe the 
electric potential and the plasma velocity fields, is solved in closed 
form. The electric field, current density, and velocity distributions 
are discussed in terms of the Hartmann number H and the Hall coeffici- 
ent (i)T. As a result of the Lorentz forces, the plasma rotates with 
speeds as high as 10^ cm/sec around its axis of symmetry at sufficiently 
large values of H and WT. It is remarkable that the Hall effect 
supports the plasma rotation. As a result of the centrifugal forces 
(in the system of reference rotating with the plasma), the heavier 
(lighter) atom and ion components are enriched in the peripheral 
(central) region of the discharge centrifuge. 


*Supported in part by NASA. 
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An electrical discharge in a cylindrical container rotates if the 
Lorentz force has a nonvanishing component in the azimuthal direction. 

For example, arc experiments in an aclal external magnetic field B 

z 

1 2 ) 

indicate that the discharge plasma rotates—’— since the current field 
lines j have a nonvanishing radial component so that (j ~ 

-j B 0. In a stable arc dlshcarge, the j -component is caused by 

IT Z "C 

the concentration of field lines (j) at the electrodes and a dilatation 

(repulsion of currents in the same direction) of the field lines j in 

the Interelectrode space.— ^ In an unstable arc discharge, significant 

radial current components j are associated with "the m = 0 and higher 

^ } 

2 ) - ‘ 

order (m ^ 1) instabilities.— 

A simple model for the production of -at?.' electrical discharge 

centrifuge, which has a radial current density which is in magnitude 

comparable with the axial current density j , is shown in Fig. 1. The 

radial spreading of the current field lines ^ is forced by means of 

electrodes of considerably different radii and (R2 >> in 

the end plates z = ±c of an electrically isolating discharge chamber 

of radius R . The field lines of the current density j and of the 
o 

external axial magnetic field cross under a nonvanishing angle 

(except at the chamber axis) so that the resultant Lorentz force j x 

rotates the discharge around its axis of symmetry. In steady state, 

the magnetic body forces in the azimuthal direction are balanced by 

the viscous forces (boundary layers at the chamber walls). 

In the literature, the observed rotation of various electrical 

discharges in axial magnetic fields appears to be understood qualitatively 

12 ) 

in simple cases.—’— Exact solutions for plasma centrifuges are apparently 

5 6 ^ 

not known either in the case of collisionless or collision-dominated— — 


plasmas . 
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In the following, the steady state rotation of the spatially diverging 

discharge contained by an insulating cylinder in the external axial magnetic 

field ^ (^'Ig' 1) is treated theoretically. The analysis is based on the 

o 

magnetogasdynamic approximation, in which two characteristic nondimens lonal 
parameters occur, the Hartman number H and the Hall coefficient ujt, 

H * (o/u)^B R , WT •* (lelB /m)x. 

The symbols designate the electrical conductivity (a), the viscosity (y), 
the electron gyration frequency (ta) , and the electron momentum relaxation 
time (t). H and wx are a measure for the strength of the Lorentz force 
relative to the viscous force and for the reduction of the current flow 
transverse to the magnetic field B^, respectively. The magnetic fields 
associated with the discharge currents * {j , j", j } are neglected for 

IT 0 Z 
2 2 

small magnetic Reynolds numbers [o_^ ■ o/(l + w x )], 

- 0[B^/B^] «» y^tDToj_vr'<<l , 

Rq = 0[Bg/B^3 ' |y^l/2TTR^l/B^«l . (1) 

R * 0[B /B ] * y wxaivR <<1 ■ , 

z z o • o o 

where v is the characteristic velocity of rotation and I the discharge 

current. These inequalities are satisfied in many cases, e.g., if 

R„«l for i) o)T>>l and R<1 or ii) R«1 and 0 < tor < «>, where R = y ov max 
0 o 


PTTPRODUCIBILITY of tu- 
WAL PAGE IS POOE 
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THEORETICAL FORMULATION 


For a purely azimuthal flow field, v ■ {0, v(r,z),0}, the plasma 
behaves incompressible, V * v * 0. From the continuity equation, 

V . (pv) = V • 7p « 0, it follows then that the density gradient Vp is 
everywhere perpendicular to the flow field v. These ideal conditions 
are realized if secondary flows are absent or at least negligible.—^ 

ON 

In accordance with the magnetogasdynamic equations,— Ohm’s law with 

Hall effect,—^ and the conservation equation for the electric charge 

density (V • ^ ■ 0), the rotating discharge in a homogeneous magnetic 

field ^ (Fig. 1) is described by the following boundary-value problem 
0 

for the azimuthal velocity v(r,z) and electric potential <(>(r,z) fields 
[induced magnetic fields neglected, Eq. (1)]: 




r 3r 


Oj_ 3z 


( 2 ) 

(3) 


where 


v(r,z) _ *=0,-c<z< + c 

’ r“R ' _ _ 


(4) 

0 



v(r,z)z.±c “ 0 , 0 < r < R^ 

9 

(5) 


and 

- a[3<|)(r,z)/3z]^^^^ = I 6(r - R^ ’ 0 1 ^ , (6) 

[3<[)(r,z)/3r]^_j^ = 0 , -c £ z _< + c . (7) 

o 

The boundary conditions (4), (5), and (7) consider that the plasma 
does not slip at the walls r = R and z = ±c, and that no current 
flows into the cylinder wall r = R^, respectively. The boundary con- 
ditions in Eq. (6) imply that the cathode (R^^) and anode (R 2 ) are ring 
electrodes of vanishing radial width, Ar ->■ 0[6(r - R^ = radial 

Dirac function] . The net current flowing through the discharge is by 


Eq. C6) 
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(5(r-R^ 2 ) 


I < 0 , 


since the positive current (I < 0) flows from the anode to the cathode 
(Fig. 1). The pressure distribution p * pCr,z) is determined by the 


r- and z-components of the equation of motion, 

■"m ^ ' I? "-L ®o<- If + ®o> 


( 8 r) 


According to Eq. ( 82 ), it is 3p/9z » 0 if the induced fields B and B- 

r D 

are neglected [Eq. (1)]. This means that momentum cannot be exactly 

balanced in the axial direction if induced magnetic fields are neglected 

(in absence of secondary flows).— ^ Eq. ( 8 z) is in accord with the boundary- 

layer approximation^^ according to which the normal pressure gradient is 

3p/9z * 0 at the electrode plates z ■ ±c, 0 < r < R . 

— — o 

In absence of the Hall effect, ojt<< 1, it is Vx^ « ^d^^r’ '^z^* 

Hence, B •* 0 and B “0 because of the homogeneity of the boimdary con- 

IT 2 

ditions for B and B , whereas 7 ^ 0 since 2 ^ 0[B„(r,z = ±c) = 

r z 0 ^ 2 o 

(p l/2ur)H(r-R_ t)]. Consideration of the induced field B = {0,B„,0} 
o z , i. 0 

leaves Eqs. (2)-(3) unchanged. This means that the boundary value prob- 
lem in Eqs. (2)-(7) and the solutions v(r,z) and ^ir,z) derived from it 
remain valid even in presence of a significant induced field B = {0,B.,0}, 

D 

R->1, as long as the Hall effect is negligible, wt«1. 




ANALYTICAL. SOLUTION 


The characteristic nondlmenslonal parameters of the magnetogasdynamic 
discharge problem under consideration are obtained by Introducing the 
nondlmenslonal Independent and dependent variables, 


p “ r/R , 0 < p < 1 

o _ — 

C « z/c , -1 £ C <. +1 


and 


where 


V(p,?) « v(r,z)/v^, $(p,5) » 4>(r,z)/(j)^ ^ 


V = 4i /R B , (|> = Ic/2iTaR^ 

O ^0000 o 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


In terms of the nondlmenslonal space variables and fields, the 
boundary-value problem defined in Eqs. (2) - (7) assumes for V(p,^) 
and $(p,C) the form: 


1 ( 3£ 

p 3p 3p 


(p_)+^2_= (pV) 


p Bp 


where 


and 


^[i^(pv)l +n-2 0-h|v = -h|_|1. 


° » -1 < ? < +1 
P-1 — - 

V(p,C)^^+3^ “ 0 , 0 < p < 1 


-[3$(p,c;)/3C3^_^^ = 5(p-P2,l)/P , O^p^l 

t34»(p,;)/3p]p^^ =0 , -1 £ C 1 1 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


with P2jl = nondlmenslonal constants H, N, and Hj^ are 

defined by 
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( 19 ) 


■ (1 + u^t2)(R^/c)^ , N“^ « (R^/c)^ , 

H]_.(a 

In view of the similarity of the left sides of Eqs. (13) - (14) 

with Bessel's differential equation, z" + p"^Z' + (k^ - p"^m^)Z ■ 0, 

ID in V ID 

for cylinder functions ^jj^(^yp)» partial solutions of the coupled 
inhomogeneous equations are sought in the form, 

4>^(p.;) - J^(k^p)f^(C)~^ (21) 

V^(P.;) » J.i(k^p)g^(!;) , (22) 

where J^(k^p) - -Ji(k^p) and J^(k^p) + (k^p)"^ Jj (k^p) = J^(k^p). 
Substitution of Eqs. (21) - (22) into Eqs. (13) - (14) yields 

d^f 

m 2 - k^ m2 . C23) 

- (k2 + Hpn2g^ . hJn2£^ . (24) 

where the eigen-values k^ > 0 are determined by the boundary conditions 
(15) and (18) as the real roots of the transcendental equation, 

Jl (k^) * 0 , V = 1,2,3,.... (25) 

Thus, the general solution of the coupled equations (13) - (14) obtains 


by linear superposition as: 


00 

$(p,c) = -2C + E J^(k^p)f^(c) , 

(26)' 

V=1 


00 

V(p,c) = B J^(k^p)g^(C) . 

(27) 


v=l 


In view of Eq. (25), Eq. (26) is a Eourier-Dini series in which a zero- 
order term, -2c, has to be included, in accordance with the Fourier-Dini 
expansion [Eq. (32)] of the boundary value in Eq. (17), whereas Eq. (27) is 
a Fourier-Bessel series. By decoupling Eqs. (23) - (24) one finds for 
f^(C) and g^(?) the differential equations of 4th order, 


)/p)b2r2 - h 2/(1 + w^T^) . (20) 

J_ o o 
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(28) 




d^f 


dC 

d^g, 


IT - + N^) + N^H^l +1 cVn^ » 0, 


d^8^ 422 


^ - [kJCM^ + n 2) + N^H^] + k^M N g^ - 0, 


(29) 


with 


- [df^(C)/dC]^.^^ - • 


(30) 




(31) 


as boundary conditions by Eqs. (16) - (17). In deriving Eq. (30), the 
Dirac function in Eq. (17) has been expanded»as the Eourier-Dini series, 

00 

^(P “ P2,i)/P“2 '+ 2 I [J (k P2 i)/J^(k )]J (k p) . (32) 

’ ,ov» ov ov 

V“1 

In addition to Eqs (28) - (31), f^(C) and g^(?i) have to satisfy also 
the uncoupled Eqs. (23) - (24). With 


0), = 0) , , w. = w , (L' = “W ,, , 0), = “W ; 

Iv V+ 2v V“ 3V V+ 4v V" 


(33) 


+ N^) + ± {[k2(M^ + N^) + 

- , (34) 

the general solutions for f^(5) ~ and g^(4) ~ e^^ of Eqs. (28) 

(29), can be written as: 


sinh w. 4 cosh io, Z 

f (C) = A. — : - ^y - + B. r— ^ 

\) Iv sinh 05, , Iv cosh 05 


Iv 


Iv 


sxnh u)„ t cosh u)„ L 

, . 2v^ . „ 2v^ 

+ A. ■ — r“T — + B 


2v sinh 05 


2v 


2v cosh 05„ ’ 

2V 


(35) 


sinh 05, 5 cosh 05, ? 

g (5) = c — r-r — ^ + D — — ^ 

°v IV sxnh 05, . IV cosh 05 


Iv 


iv 


sinh 05„ c cosh 05_ t 

+ C. — — - + D 


2v sinh 05^^ 2v cosh 05^^ 


(36) 
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Only four of the eight integration constants for any 

V ^ 1 are independent; by Eqs (23) - (2A) , 

(u)f - k2M2)A, - k M^C, 

' Iv V Iv V Iv 


(to; , - k2M2)A^ - k 

2V V ' 2V V 2V 


(to? - k2M2)B, - k M^D, 

' Iv V Iv V Iv 


(37) 


(toj - k2M2)B- ^ - k M^D„ ^ 
2V V 2\) V 2\> 


and 


h|^)n2]c 


iv 


[to^ 

^ 2V 


t“lv 


[“iv 


(k2 + h[)n2)C^^ 
(kj + h[)n2ID;^ 
(kj + 


k N^H?A, 

V _L Iv 

k N^H? A, 

V J_ 2v 

k N^H?B, 

V I v 

k n2r?B^ 

V J_ 2v 


(38) 


where the coefficient determinants of the pairs of corresponding 
equations in Eqs. (37) and (38) vanish owing to Eqs. (33) “ (34). 

Upon application of the four relations in Eq. (38) , which are 
equivalent to Eq. (37) by Eqs. (33) - (34), and the boundary conditions 
(31) , which give 


“ c C , - c 

2v Iv V 


-D 


2v Iv 


= D 


(39) 


Equations (35) - (36) become: 


' k N^H? j'^lv slnh u 

V i ( 


sinh 0), C sinh (o„ t 

IV^ . „ 


Iv 


2v sinh to 


2v 


j ^ cosh to^^C cosh 

^ k^N^H|^ I ^Iv cosh to^^ ^2v cosh to^^ 


(40) 
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-^v 


g (c) ■ c 

®V V 


r sinh sinh 


slnh 0) 


iv 


sinh 

sinh 


+ D 


fcosh cosh w 


cosh u 


Iv 


cosh 


kail 

“ 2 vJ ■ 


( 41 ) 


where 


“iv = “iv - <“5 i - 1-2 • 

The boundary conditions (30) applied to Eq. (40) yield 

c ■ Jq<V 2 >’ 

'“iv^lv'^^lv -“2 v"2v“*'“2vJ 


(42) 


(43) 


k N^H| 

1 


~ -^0 ^2^^ 


,t“lvSv^8h«-lV " 


(44) 


Substitution of Eqs. (43) - (44) into Eqs. (40) - (41) gives as 
solutions for f^(?) 3^^ g^(?)! 




‘ ■ [u)^^f2j^ctho),.. - w,. cthu^..] 


Iv 2V 2V 2V 




sinh a), C 
Iv 


Iv sinh 03 


- n 


sinh 03^^? 


iv 


2v sinh 03 


2 v J 


l“lV^lv'«*“’lv'“2v"2v‘«'“2v' 


x[(!l 


cosh 03, 5 
Iv 


V cosh 03 


- 0, 


iv 


cosh 
2 V cosh 


“2V-I ’ 


(45) 


and 
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11 


) 

gXO • 7-7 




+ ^0<V2» 


■sinh sinh 

slnh u... ” sinh w 


IV 


2v J 


t-'o<Vl> - Jo<V2>i 


"cosh cosh 

I cosh 0 ), 


cosh u). 


(46) 


'iv — "av -• 

Equations (45) -- (46) form, together with Eqs. (26) - (27), the closed 
form solution of the problem of the rotating gas discharge in an axial 


magnetic field B^: 




* J (k p) 

-tfo C) - - 2c - 7 0 

iP.u I I ^ n cthu), -oj^ cthw- ] 

V"1 j"^(k ) ( '■ Iv Iv iv 2V 2V 2v" 




O V 


sinh sinh "| 


V sinh u 


2 V sinh 0 ), 


IV ■ 2v 

tJptVl) - Jq<V 2>^ 


Q 


cosh 0 ), ? 
Iv 


Iv cosh oj 


n 


IV 


cosh ■ 

2 v cosh "w^” 

(47) 


and 


; Vi<V> 


O V 




v(p.c) - - j 2 (jj % I [u,..!l,.,cthu,„ -u.„!l,, cthu.„J 


'Iv iv 


1 V 2V 2 V 


fsinh w, ? sinh C 
' iv 2v 

sinh 0 ), . sinh w 


Iv 


2V -J 


[J (k p ) - J (k p,)] 

'■ o V 1 o V 2 


’cosh cosh 

cosh M. 


IV 


cosh ( 1 ) 


2V -I 


(A8) 


["lv"lv'^8ho.,„ - 032^n2^tgto2j 

The remaining nondimensional discharge fields E = - 

? = J^/J’ are given in terms of the solutions for $(p,C) and V(p,?): 


- 1 , 


E = a$/3p , E. = 0 , E = -N ’^3$/3S 


(49) 
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m 


12 


J 

r 




1 + 






COT 


1 + 


CO^T' 




J = 
z 


1 ^ 

N 

(50) 


where E = <}) </R » j * cr<<> /R * and N “ c/R [Eq. (12)]. 
o o o o o o 

If the cathode is In the plane z = -c(C = ”1) and the anode is 
in the plane z = +c(? « +1), then the reference fields and 

[Eq. (12)] are negative, since I <0. The results are also applicable 
to the case where the anode is in the plane z = -c(C = -1) and the 
cathode is in the plane z = +c(c; = +1). In the latter situation, the 
reference fields and <f>^ [Eq. (12)] are positive, since I > 0. 

These explanations hold for magnetic fields pointing in the positive 
z- direction, B > 0; v changes its sign if B < 0 [Eq. (12)]. 


..^EODOCIBILmOF™ 

I ’.ifNAL PA(jE IS POOR 


^48 


13 


NUMERICAL ILLUSTRATIONS 

As an Illustration the radial (p) dependence of the nondlmenslonal 
discharge fields V(p»C), E^XPif.), E^(p,C)» «nd J|,(p,c) luis 

been calculated for I <0 In the cross-sectional planes C = -0.99 
(cathode region), C = 0 (central region), and ^ = +0.99 (anode region) 

based on Eqs. (47)-(50). The remaining fields J„(p,?) and J (p,?) 

t) z 

are proportional to J (p,?) and E (p,C), respectively [Eq. (50)]. 

r z 

The characteristic (nondlmenslonal) magnetic interaction numbers are 
treated as parameters: 

OJT “ 1, 10 ; H = 1, 10, 100. 

-2 2 2 

The geometry parameter N is taken to be N = 1 so that M = 1 + w t , 

corresponding to R^ c [Eq. (20)]. The radial positions of the 

cathode and anode are assumed to be: 

p, = 0.01 (R, = 0.01 R ) ; p- = 0.9 (R_ = 0.9 R ). 

"^1 1 o ' ^ ^2 2 o 

The dimensional fields are negative everywhere where the nondlmenslonal 
fields are positive, and vice-versa [Eq. (11) J since < 0 and 
< 0 for I < 0 [Eq. (12)]. 

Central Region, ^ = 0 : In the Figs. 2-6, the azimuthal velocity field 

V(p,0), the electric potential $(p,0), the radial and axial electric 

fields E (p,0) and E (p,0)«J (p,0), and the radial current density 
r z z 

J^(p,0) are represented versus p, with (ux,H) = (1,1), (1,10) , (1,100) , 
(10,1), (10,10), (10,100) as parameters. It is seen that jv] 
increases considerably at any point 0 < p < 1 if either H or wt 
are increased. In the region p > 0 sufficiently close to the axis, 

J$j, |e 1, ]e - 2N ^1, and jj | increase with increasing H or wx. 

The field distributions move towards the axis p = 0 as cox becomes 
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larger. The "hump" developing at p = 0.9 (Figs. 4-6) with increasing 
(ijT shows the influence of the ring anode (p = 0.9, C = +1) in the 
plane i; " 0. 

Cathode Region, C » -0,99 ; The Figs. 7-11 show V(p, -0.99), 

$(p, -0.99), E (p, -0.99), E (p, -0.99)“J (p, -0.99), and J (p, -0.99) 
versus p with (ujt,H) « (1,1) ,... (10,100) as parameters. These fields 
increase in intensity at any point 0 < p < 1 if H or wt is 
increased. Since the ring cathode is at p^^ = 0.01 (? = -1), the field 
distributions are closer concentrated at the axis p = 0 than those 
in the plane 0 (Figs. 2-6). Note that the plasma rotates only 

in the region p = 0.1 with a significant velocity, since the 
Lorentz force decreases rapidly with increasing p-^l. A 

comparison of the corresponding fields in Figs. 2-6 and Figs. 7-11 
indicates that the discharge spreads slightly in radial direction 
with increasing -1 < ? £ 0. In particular, an increasing radial 
section of the plasma rotates with a significant speed as -1 < ^ ^ 0 
increases. 

Anode Region. C ~ +0.99 : In the Figs. 12-16, V(p, +0.99), ^(P> +0.99), 
E (p, +0.99), E (p, +0.99)«J (p, +0.99), and J (p, +0,99) 

IT Z Z ^ 

are plotted versus p with (o3t,H) = (1,1) ,... (10,100) 
as parameters. The dependence of these fields on H and wt is as 
in the previous cases for ? = 0 and % = -0,99. The velocity 
distributions are fully developed nearly through the entire chamber 
across section 0 < p <0.9, since the Lorentz force -J^B is 
strongest in the vicinity p = 0.9 of the ring anode, p = 0.9(C = +1) . 
As a result, a thin and steep boundary layer exists close to the 
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cylinder wall (p * 1) with backflows at sufficiently small wx-values 
(Fig. 12) ^ The radial distributions of E , E “:J , and J (Figs. 13-16) 
clenrly Indicate that, in the plane C " +0.99, the electrical dlachorge 
has shifted to the region p = 0.9 due to the influence of the 
nearby ring anode, p ■ 0.9(C * +1) • This shift occurs first slowly in 
the region -1 < 5 < +1 - A?, and then rapidly in a relatively thin 
layer << 1 close to the anode plane ? = +1. 

It is remarkable that the discharge remains concentrated in a 
radial region close to the cylinder axis with little radial spreading 
of the current density except in a layer AC close to the ring 
electrode of large radius (R 2 » R^) in which the radial current 
component dominates the axial current component J^. This 

spatial concentration of the discharge 1s the more pronounced the 
larger H and 0 )t are, since the axial magnetic field reduces 

the radial current flow for increasing ux. The speed of 

plasma rotation V(p,C) increases with increasing magnetic induction 
by orders of magnitude over the reference speed as the 

Figs. 2, 7, and 12 indicate which show V(p,C) for increasing wx 
and H. The theoretical electric field and current density distribu- 
tions are in qualitative agreement with experiments.—^ 

The graphs in Figs. 2-16 are based on the Fourier -series solutions, 
in which the first 100 terms were considered and the eigenvalues 
were calculated up to the 10th decimal point. An even larger number of 
terms in the Fourier series solutions has to be taken into account if 
one wishes to compute v pproximately) the discharge fields extremely 
close to the ring cathode (p = 0.01, C ^1) and ring anode 
(p = 0.9, ? = +1) where 9$(p,?)/95 changes discontinuously with p 
due to the electrode boundary conditions. 


- 51 - 


16 


APPLICATION 

The syetero analysis presented indicates that extremely high 
speeds of plasma rotation are obtainable already at moderate discharge 
currents I and magnetic inductions presumed the magnetic inter- 

action numbers are not small, H > 1, ujt > 1. As an example, consider 
a rotating arc discharge with: 

|lj = 10^ amp , |b^| = 10° Tesla, 

2 -1 -1 
a = 10 mho m , R = c = 10 m . 

* o 

Hence, by Eq. (12) 

3 1 —1 

V " Ic/2iraB R " (5/ir) x 10 m sec , 

0 o o 

and, by Fig. 2, 

0[v] ■ 0[v^V] » 10^ m sec for wt = 10, H = 100. 

Speeds of plasma rotation v, which are by orders of magnitude larger 
6 ""1 

than 10 cm sec , can be produced if the order of magnitude of the 
parameters wt and H is increased. The viscous forces reduce, 
however, the speed or plasma rotation always in the layers close to 
the walls (z = ±c; r = ' 

One recognizes that the torque of the Lorentz forces is a large 
physical effect which can be readily used to build an electrical 
"discharge-ultra-centrifuge" for isotope separation. If the working 
gas of the electrical discharge consists of two isotope gases, then 
the centrifugal forces would concentrate the lighter isotope ions and 
atoms in the central region and enrich the heavier isotope atoms and 
ions in the peripheral region of the discharge. 


- 52 - 


17 


As an illustration to the plasma centrifuge, the isotope 
separation ratio Is calculated. According to the equations of motion 
for two isotopes of masses and m^ , the Isotope density ratio at 

distances 0 £ r £ ~ where 6 is the viscous boundary layer 

thickness, is approximately (T^ = temperature of the isotope ions) 

n, (r) ^ n. (0) -»^Am .v(r)^/kT 

= e ^ , Am = m . - m , 

Hj(r) n^iO) 13 X :) 


where the bar designates a spatial average over the region Iz] < c. As 

a specific example, consider an uranium plasma centrifuge containing the 

isotope ions (1) and (j) at a temperature = 10 "K 

(and electrons at a temperature > T^) • In this case, one has 

-27 -20 

Am^j * m(237) - m(235) = 3.320 x 10 ^'kg, kT = 1.381 x 10 Joule. 
Hence, the Isotope separation ratio Is: 


n„«.,(r) n„o-^(0) __ n _ 3 _i 

/ _£dL 1.X28 X 10 for v(r) = 10 m sec , 

n23s(r) ’^235^®V 

=1.661 X 10^ for v(r) = 10^ m sec 


It should be noted that the magnetogasdynamic considerations are 

only applicable to dense discharge plasmas, for which the mean free 

paths of the electrons, ions, and atoms are small compared to the 

characteristic chamber dimension, min (R^;c). The evaluation of collision 

less, rotating plasmas poses a rather different problem, since in this 

case the interactions through the self consistent electric magnetic 

9) 

fields play a dominant role.— The separation of isotopes by centrifugal 
forces in low density plasmas has been established experimentally.—^ 
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fig. 5: E^(p,c) versus p for S = 0, and (wt, H) = (1,1) to (10,100) 
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J (p,C) versus p for ? = -0.99, and (ojt, H) - (1,1) to (10,100) 











FIG. 16: J^(p,i;) versus p for ^ - +0.99, and (tux, H) = (1,1) to (10,100). 
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DEPOSITION OF SPUTTERING PRODUCTS ON SYSTEM SURFACES* 


H. E. Wilhelm and S. H. Hong 
Department of Electrical Engineering 
Colorado State University 
Fort Collins, Colorado 80523 


ABSTRACT 

The boundary-value problem describing the diffusion of sputtered 
atoms and their deposition on the surfaces of ion propulsion systems is 
treated by means of a model of simple geometry. Based on an integral 
theorem of Weber for outer boundary-value problems with an inner 
cylindrical boundary, an analytical solution is derived for the density 
field of the sputtered atoms which involves a function which is deter- 
mined by an inhomogeneous integral equation. 


^Supported by NASA 
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In an ideal vacuum, sputtered atoms travel undeflected along 
straight paths determined by their initial velocities at the point of 
emission. Within this free particle flow, a system surface is reached 
by the sputtered atoms only if it can be seen along a straight line from 
the emitting surface. In reality, ion propulsion systems are surrounded 
by a very rarefied plasma consisting of escaped beam ions, recombined 
ions, and electrons. For this reason, always some of the sputtered 
atoms will be deflected out of their initial paths by interacting through 
long-range forces (polarization forces) with the plasma particles so that 
they can reach system surfaces which are not seen along a straight line 
from the emitter. 

Figure 1 depicts the geometry of an idealized propulsion system 
which exhibits an emitting p.lane z»0, O^r^a (accelerating grid), 
the rocket surfaces r«a, -c<^^j<0 and z = -c, 0 1 r ^ a, and the 
plane z = -d, a t _< b of the solar energy collectors. All these 
system surfaces can be reached by the atoms sputtered from the emitter 
by diffusion through the rarefied plasma. The diffusion coefficient 
D is determined by the Vlasov equation^^ for the sputtered atoms inter- 
acting through weak long-range forces^^ with the plasma particles . In 
view of the mathematical difficulties associated with the solution of 
outer boundary-value problems for the geometry in Fig. 1, a somewhat 
simpler system is studied here consisting of an emitting plane 
(z = 0, 0 £ r ^ a) , the upper rocket surface (r = a, -c £ z £ 0) and 
the plane (z = -c, a £ r £ >») of the solar energy collectors (Fig. 2). 

The latter is assumed to have infinite radial extension, r = b , 

max 

since in general b>>a and b>>c,d (Fig. 1) . Within the model of 
Fig. 2, particle deposition on system surfaces in the space z £ -c 
cannot be analyzed. 
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BOUNDARY-VALUE PROBLEM 


In the space z -c, let the density of the sputtered atoms be 

_3 

designated by n(r,z) [cm ] and the flux of emitted atoms at the 

-3 -1 

emitter surface by I(r) [cm • cm sec ]. In steady state, the 
spatial distribution n = n(r,z) of sputtered atoms is determined by 
the boundary-value problem for the Laplace diffusion equation (Fig. 2) 


where 


3^n . 1 ^ ^ 


[9n(r,z)/9z]^^Q = -I(r)D 


0 < r < a 


and 


n(r,z) =0, -c < z < 0 
» 'r=a * _ _ 




2 2 

n(r,z)->0, (r +z)->«> 


( 1 ) 

( 2 ) 

, (3) 

, (4) 

, (5) 


are the proper and improper boundary conditions, respectively. D 
designates the diffusion coefficient of the sputtered atoms in the 
rarefied plasma which represents a spatial average, D = < D(r,z) >. 

The boundary conditions (3)- (4) imply that sputtered atoms 
arriving at the system surfaces are deposited there, i.e, , do not 
return into the diffusion space. This assumption is at least 
approximately correct for nonheated surfaces as long as the number of 
atomic layers deposited is not too large. The fluxes = -D V^n 
of atoms arriving at the system surfaces r = a, -c ^ z _< 0 
z - -c, a ^ r ^ are given by (Fig. 2); 


and 
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41 ^ (r • a,z) •“ -D On(r “ a,z)/3r, -c £ z £ 0 , (6) 

4> (z » -c,r) » -D 3n(z «» -c,r)/3z, a _< r £ »> , (7) 

Accordingly, 


0 

_ = -2iTaD / [3n(r = a,z)/Srldz , (8) 

-c 

00 

^z=-c ~ “27 tD J [3n(r,z = -c)/3z]rdr , (9) 

a 

are the numbers of sputtered atoms deposited per unit time on the 
system surfaces r = a, -c £ z £ 0 and z = -c, a ^ r ^ respectively. 
The above boundary-value problem cannot be solved directly, i.e., 
requires a decomposition of the space z ^ “C into appropriate 
subregions for which the associated boundary-value problems are 
solvable. In this approach, the common boundary value [4^(r)] at the 
decomposition plane is detexrmined by an integral equation. 

Let nondimens ional independent and dependent variables 
be introduced in accordance with: 

p = r/a, 0 _< p £ C = z/c, "1 ^ ^ “ » (1^) 

and 

N(p,5) = n(r,z)/n^, S(p) = I(r)/l^ , (11) 

with 

n^ = cI^/D, = l(r=0), y = c/a . (12) 

The boundary-value problem defined in Eqs. (l)-(5) reads in 
nondimensional form: 


9 N 1 


3N 


3p 


P 9p 


+ Y 


2 32n 




(13) 
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where 


[3N(p,C)/3?]^^q « -S(p), 0<p<l 


N(p,5)p,^ “0, "1 1 5 1 0 


N(p..;)^^_j^ - 0, 1 5 p 1 


N(p.?) 0, (p^ + 


In Fig. 2, the space is decomposed into the regions 
1(0 £ p jS «®» 0 ^ and 11(1 ^ P ^ "j “1 1. C ^ 

interface ^ = 0, 1 ^ p ^ the partial 3N(p,^ = 0)/3C = i{;(p)H(p-l) 
is introduced as the common (unknown) boundary value i|^(p) of the 
adjacent regions I and IX, 1 5. P ^ '^^ns, the boundary-value problem 

in Eqs. (13) -(17) can be decomposed into boundary-value problems 
for the regions I and II. 

I. In region I, N = Nj(p,?) is described by the ’’regular" 
boundary-value problem: 


, 1 ^ -2 a^N 


0, 0<P<“>, 0<c< 


[3Nj(P,0/3C]^^q = -S(p), 0 < p < 1 - 0, 

= Y(p) , 1 + 0 £ p 5 “ 


where 


N^(p,?> 0, (p^ + C^) -> 


S(p=l) = 4'(p=l) = 0 


, ( 20 ) 
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for physical reasons. Since region I is the upper half of the 
infinite space (0 £ C £ «) , the general solution of Eqs. (18) and 
(20) is given by the Fourier integral, 

00 

N (p,C) ■ f A J (kp)dk , (22) 

0 ° 

which satisfies the improper boundary condition for p ^ and 
C The Fourier amplitude A(k) is determined by the boundary 

condition (19) , 


00 

-Y / A(k) J (kp) kdk = - S(p)H(l-p) + 'i'(p)H(p-l) . (23) 

0 ° 

Application of the inverse Hankel transform to Eq. (23) gives; 

- 1 

A(k) = y"^ / S(a)J (ka)ada - y”^ / ^(ct)J (ka)ada , (24) 

0 ° 1 ° 

Substitution of Eq.(24) into Eq. (22) results In the solution for 
region I; 


- 1 1 ^ 

N-(p,C) = Y / 'ike ^ ^ J (kp)[/ S(a)J (ka)ada -- / ^i^(a)J (ka)ada] > 
^ 0 ° 0 ° 1 ° 


0<p<«, 0<?< 


(25) 


II. In region II, N h N^j.(p,^) is described by the "outer" 
boundary-value problem: 


^ ^TI 1 ^^II -2 ^ ^11 

^ + ^=0, 0 < p < CO, -1 < c < 0 

3p^ P 3p 3^2 

[3Nj^(p,?)/3?]^^^ = 'i'(p), l£P£“ 




Nii(p,?)p^l= 0, -K5<0 


(26) 

(27) 

(28) 




-78 



9 


Nii(p,C)^._l - 0, 1 1 P < ® » (29) 

Njj(p.C) 0, P -1 £ ; <. 0 . (30) 


According to Eq. (28), region II has an inner, cylindrical boundary 


at p = 1 where N^j(p,?) vanishes. For this reason, a Fourier 

integral representation of Nj.j(p,C) is needed for 1 £ p ^ ® 

which vanishes at p “ 1. According to Weber's integral theorem, an 

arbitrary function ij) (p) , <2 £ p £ with \|;(p = a, “) = 0 satisfies 

2 3) 

the integral equation: — 


'i'<p) 


kdkW^(p) 


Q jJ(ka) + Y^(ka) 


ijj (a)Wj^(a)ada, a £ p £ “ 


(31) 


a 


where 


Wjj(p) S J^(kp)Y^(ka) - J^(ka)Y^(kp) , (32) 

and J^(kp) and Y^(kp) are Bessel functions of order v of the 
first and second kind, resepctively . In view of Eqs. (31)-(32), a 
Fourier integral solution of Eqs. (26)-(30), is sought in the form, 

00 

Nj^(p,C) = / B(k)W^(p) slnkYk(i;+l) dk , (33) 

W^(p) 5 Jo(kp)Yo(k) - Jo(k)Yo(kp) , (34) 

which obviously satisfies the boundary conditions (28)- (30). The 
Fourier amplitude B(k) is determined by the boundary condition (27), 

, (35) 

, (36) 


CO 

Y / B(k)W. (p) coshYk kdk = H'(p), 1 £ p £ “ 

0 

which gives 

00 

B(k) = Y ^cosh ^Yk/ H'(“)W^^(tt)oida/[J^(k) -f Y^(k)] 


-79 


10 


by Eq. (31). Substitution of Eq. (36) into Eq. (33) results in the 
solution for region II: 


N^.j(p,C) 



i coshyh 


1 <_ p 


W, (p) " 

“5 7— J H'(a)W, (a)ada, 

r(k)+Y^(k) 1 

-1 < ? < 0 


(37) 


The solutions 'N^(p,C), Eq. (25) , and N^^(p,^) Eq. (37), covitain 
the yet unknown boundary-value 'l'(oi), 1 ^ a ^ 'l'(a) is determined 

by the continuity condition at the interface of regions I and II, 


Nj (p,e)^^Q =“ Nj^(p,?)^^q, 1 1 P 1 “» > (38) 

which gives 

00 2 . 00 

JdkJ (kp)[Js(a)J (ka)ctda - J'i'(a)J (ka)ada] 

0 ° 0 ° 1 ° 

= Jdktghyk — r x /4'(a)W, (a)ada, 1 £ P ^ “ • (39) 

0 J^(k)+Y (k) 1 

o o' ' 

Eq. (39) indicates that 'i'(a) is determined by an inhomogeneous integral 
equation.—^ 

Because of the boundary conditions (19) and (27), the remaining 
continuity condition at the interface of regions I and II, 

[3N^(p,?)/a?]^^Q = [3 N^j(P,?)/3C]^^q, 1 1 P 1 “ , (AO) 

has already been satisfied. Indeed, substitution of Eqs. (25) and (37) 
into Eq. (40) yields 
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-/ kdkJ (kp)[/S(a)J (ka)ada - Jf(a)J (ka)ada] 

1 ° 


0 ° 0 ° 


/ kdk -= r / H'(a)W, (a)ada, 

0 J'^(k)+Y (k) 1 

o o 


1 < p < 


Since, 


Jj (kp)J (ka) kdk 
0 ^ ° 


= 6(a-p)/a 


~ W^(p)Wj^(a) 

/ 5 kdk 

0 r(k)+Y^(k) 

o o 


= 6(a-p)/a 


Equation (Al) reduces to the expression [6(a-p) = Dirac function] 


- /s(a)6(a-p) da + /'l'(a)6 (a-p)da 
0 1 


= /H'(a)6(a-p)da, 


1 < p < 


which gives the expected identity, ^ ip) = H'(p) 

INTEGRAL EQUATION 

By introducing the kernel fC(a,p) and the source Q(p) , the 
Integral equation in Eq. (39) can be rewritten in the convenient 


J'J'(a)K(a,p)ada = Q(p), 1 jl P £ 

1 


where 


“ W (a)W, (p) 

K(a,p) = / [J (ka)J (kp) + tghyk 


J^(k)+Y^(k) 
o o 


-] dk 


Q(p) " JdkJ (kp)/s(a)J (ka)ada 
0 ° 0 ° 
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(Al) 

(A2) 

(A3) 

(AA) 

form: 

(A5) 

(A6) 

(A7) 
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For simple particle emission distributions S(a), e.g., in the case of 
a homogeneous and a parabolic emission distributions, respectively, the 
source integral Q(p) is readily evaluated, 

Q(p) = I p[E(^) - (1 - ^)K(^)1, 1 < p < ~ 

P 

for ,S(a) =1, 0 a ^ 1 , (A8) 

and 

Q(p) = I piuh - (1 - - ^[(1 + 4p^)e(^) - (3 + 4p^)(i - ^)uh]} 

TT p piS p , p p ^ 

for S(a)= 1-a^, » (49) 

11 

where K(~) and E(— ) are the complete elliptic integraOis of the first and 
P P 

second kind, respectively. 

Eq. (45) reduces the deposition problem to the resolution of an 
integral equation for the unknown boundary-value *F(p), 1 ^ P ^ Once 
'Kp) is determined, the particle distributions N^(p,C) and Nj.j(p,C) are 
given by Eqs. (25) and (37), respectively. From the solutions for N^(p,?) 
and N^^(p,?) follow then the deposition rates at the system surfaces ixi 
accordance with Eqs. (8) - (9). 

Unfortunately, the resolution of the integral equation leads to 
considerable mathematical difficulties. Both analytical and numerical 
methods for solving Eq. (45) have not produced final results because 
of convergence and uniqueness problems, which are typical for integral 
equations of this type. » ' 
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